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ABSTRACT 

We present results from a spectroscopic Spitzer and Herschel mid-to-far- 
infrared study of the circumbinary disk orbiting the evolved (age ~12-23 Myr) 
close binary T Tauri system V4046 Sgr. Spitzer IRS spectra show emission lines 
of [Ne il], H 2 S(l), C0 2 and HCN, while Herschel PACS and SPIRE spectra 
reveal emission from [O i], OH, and tentative detections of H 2 0 and high-J tran¬ 
sitions of CO. We measure [Ne IIl]/[Ne II] < 0.13, which is comparable to other 
X-ray/EUV luminous T Tauri stars that lack jets. We use the H 2 S(l) line lu¬ 
minosity to estimate the gas mass in the relatively warm surface layers of the 
inner disk. The presence of [O i] emission suggests that CO, H 2 0, and/or OH 
is being photodissociated, and the lack of [C i] emission suggests any excess C 
may be locked up in HCN, CN and other organic molecules. Modeling of silicate 
dust grain emission features in the mid-infrared indicates that the inner disk 
is composed mainly of large (r~ 5 /irn) amorphous pyroxene and olivine grains 
(~86% by mass) with a relatively large proportion of crystalline silicates. These 
results are consistent with other lines of evidence indicating that planet building 
is ongoing in regions of the disk within ~30 AU of the central, close binary. 
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Introduction 


The relatively rare examples of nearby T Tauri stars (i.e., those within ~100 pc of Earth) 
are excellent candidates for studying planet formation processes. Their proximity allows us 
to analyze the contents and structure of circumstcllar disks in which young protoplanets are 
likely forming. Near to far-infrared spectroscopy allows us to probe the contents of these pro¬ 
toplanetary disks, as the strongest emission from such disks arises at infrared and (sub)-mm 
wavelengths. Mid-IR (5-40 /mi) spectroscopy traces both the gas and dust components of 
the inner (R <10 AU) planet forming regions of the disk. At medium to high resolution, we 
detect em ission features due to increasing disk temper ature from the disk midplanes to their 
exteriors (ICalvet et al.lll99ll; iMalbet fe BertoutlllOOll) . In most cases, such spectral features 
only probe the disk atmosphere, as the disk becomes optically thick towards the midplane. 
Emission features may also arise from optically thin region s of the disk due to dust sublima¬ 
tion or grain growth and planetesimal formation (Pontoppidan et al. 2014). Transition disks 
with an inner gap carved out by giant planet formation and/or disk photoevaporation may 
display cavity myalls with emi ssion features due to being directly exposed to stellar UV radia¬ 
tion (e.g. Cleeves et al. 2 0111 ). Far-IR (50-700 /mi) spectroscopy traces emission from colder 
gas and dust in the outer (R >10 AU) disk. At these longer wavelengths, we expect emission 
features from gas mass tracers such as CO and atomic species that trace photodissociation 
regions at disk surfaces. 

Previous Spitzer space telescope and ground based mid-IR spectroscopic studies of pre- 
main sequence (pre-M S) star/disk systems have revealed a variety of atomic and molecular 


emission features (e.g. ISalvkl 1201 ll: IPontoppidan et al.lI2010T) . Molecules such as C 2 H, HCN, 


CN, OH, H 2 0 and C 2 H 2 , whose abundances are sensitive to stellar UV and X-ray radiation, 
are used to gau ge the stre ngth of dis k gas p hotodissociation and radiative heating processes 

(e.g. 


Pascucci et al. 2013 . 2009 : Henning et al. 2010 ; Bergin et al. 20031) . H 2 Q has also bee n 


observed in the planet forming regions around many young stars (e.g. Salvk et al.l 2 0111). 


Metal ions such as 


2012; P ascucc i et al. 


Ne il l have also been detected towards tr ansitional disks (ISacco et al. 


20071; lLahuis et al.ll2007l: iGiidel et al.ll2010l) and may trace the influence 


of high energy radiation on disk photoevaporation rates. 

Far-IR spectroscopic observations of pre-MS star/disk systems via the Herschel space 
observatory help constrain the disk gas and dust ma ss, as well as the ef f ects of high energy 


radiation on outer disk heating and chemistry (e . g. iFedele et al 


Riviere-Marichalar et al.l 120131; iHoward et al.l 120121; iThi et al.l 120101 ) . Far-IR emission from 


2013; Meeus et al. 2013 


pure rotational CO lines probe temperature regimes between those of the near-IR vibra¬ 
tional CO lines and the low-J sub-mm/radio CO li nes and ca n reveal information about this 


distribution of gas throughout the disk (e.g. Bruderer et ah _2012). The presence of far-IR 
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[O i], [C i] and [C n] emission lines may indicate that UV or X-ray radiation is incident on 


the o uter disk and photodissociating CO or OH molecules (IMa/mon et ahlll988t lAresu et al. 


2012). Modeling these and other diagnostic emission lines in Herschel spectra, combined 


with models of Spitzer spectroscopic observations, can constrain total disk mass, as well as 
the composition and distribution of gas and dust throughout the disk. 


Here, we present a combined spectroscopic Sy itzer and Herschel study of the protoplan¬ 
etary disk orbiting the nearby (D ~73 pc; Torre s et al. 2 0081 ) pre-MS binary V4046 Sgr. 
V4046 Sgr is an isolated, spectr oscopic binary that is a member of the 0 P i c Moving Group, 
with an age ~12-23 Myr old ( Torres et al. 2006 : Rinks fc Jeffries 2013 : Mamaiek &; Bell 


20141) . The binary consists of n early equal mass compon ents with spectral type s K5Ve and 


K7Ve ( Stempels fc Gahm 2004 ) and total mass 1.75 M^ (Rosenfeld et al. 2012a ). The stars 
are separated by 0.04 AU (p ~ 2.4 days; Stempels & Gahm 20041) and are surrounded by 
a large (R ~370 AU) circumbinary disk of g as and dust inclined at ~ 33° from face on 
( Rosenfeld et al. 2012a : Rodriguez et al.l 2010). The central stars are a ctively accreting ma¬ 
terial (iStempels fc Gahmll2004l: iDonati et al.ll2011c lArgiroffi et al.ll2012l) from the disk, which 
h as an inner gap of radius ~ 30 AU that is possibly due to ongoing Jovian planet formation 
((Rosenfeld et al.l 120131) . The disk is al so evidently under going photoevaporation by high en¬ 
ergy radiation from the central stars (ISacco et al.l 120121). which, in conjunct ion with planet 
formation, may contribute to the inner disk clearing ({Alexander et al.l [2014. and references 
therein) 


In this paper, we present a census of atomic and molecular emission features detected in 
the Spitzer and Herschel spectra of V4046 Sgr, and we model the continuum of the Spitzer 
spectrum to determine the dust particle composition of the disk. Our Spitzer and Herschel 
observations and data reduction processes are discussed in Section 2 and in Sections 3 & 4 
we report the results of our emission line study. In Section 5 we discuss our modeling of 
the dust features in the Spitzer spectrum and we present a summary of our conclusions in 
Section 6. 


2. Observations and Data Reduction 


2.1. Spitzer IRS 


Low (R ~60-130) and high (R ~600) resolution Spitzer Space Telescope (We rner et a l 


2004 ) data for V4046 Sgr were obtained with the InfraRed Spectrograph (IRS; IH ouck et al 
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2004 ) in April 2005 (PI: Mitsuhiko Honda@), and were retrieved through the Spitzer Heritage 


Archive. Short low (SL, 5-20/xm), short high (SH, 10-20 um), and long high (LH, 20-37 /mi) 
data were reduced individually using SMART v8.2.5 fjHigdon et all 2004)- No background 
subtraction or low-level rogue pixel removal was performed on the high resolution data 
because no dedicated background images were available. All images were manually examined 
after processing and bad pixels (particularly near 6.9 and 17.6 /irn) were smoothed over. 


The two nods of each set of cleaned spectra were averaged together and the three 
modes were combined using custom IDL programs. SL data are available for the ~5-14 /irn 
range, but are used only in the 5-10 /irn range, where high resolution data are not available. 
Examination of the original spectra showed a relatively smooth continuum, with a small 
jump between the SH and LH spectra at ~20 /irn. As V4046 Sgr is unresolved by Spitzer , 
this mismatch is likely due to the source falling off the slit in the SH module. Hence, the 
SL and SH spectra were scaled by a multiplicative factor of 1.182 such that the SL and 
SH matched the LH spectra at ~20 /mi. The final resulting Spitzer spectrum (5-35 /mi) is 
shown in Figure [U It is possible that the excess flux in the LH spectrum is due to residual 
sky background. If this is the case, then the LH spectra should be scaled down by 1.182 to 
match that of the SH and SL spectra. 


2.2. Herschel PACS and SPIRE 


Herschel Space O bservatory (Pi 


Spectrometer (PACS; IPoglitsch et al. 


.bratt et ah 120101) Photodetector Array Camera and 
201qFT range scan spectra and Spectral and Photomet¬ 


ric Imaging Receiver (SPIRE; Griffin et ah 2010) spectra of V4046 Sgr were obtained as part 
of the Cycle 1 General Observer Progranj^ (PI: G. Sacco) in June and July of 2012. SPIRE 
data were obtained at medium resolution (R ~160), and PACS data were obtained in SED 
mode (R ~940-5500). Together, these observations cover the 55-670 /mi range. Examination 
of the SPIRE level 2 data products showed a discontinuity between the two bands (190-310 
/irn and 300-670 /an) , suggesting that the background was incorrectly subtracted by the 
pipe line!!. We therefore reduced both the PACS and SPIRE data from level 0 in HIPE vlO 
(lO tt 201Ql) using interactive background normalization scripts, which allow us to subtract 
the background more accurately than is done in the standard pipeline. Data from calibration 
trees 48 and 10.1 for PACS and SPIRE, respectively, were used to correct for instrumenta- 


1 program ID: 3580, AORKEY: 11197440 
2 Observation IDs: 1342231043, 1342242448, 1342242449 
3 See section 6.4 in the SPIRE data reduction Guide v2.1. 
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tion effects and to flux calibrate the spectra. The center spatial pixel (spaxel) was extracted 
from the resulting 5x5 spaxel array and order edges and regions in the spectrum where light 
leakage occurs were removed. The resulting Herschel PACS and SPIRE data are shown in 
Figures [2] and El 

Herschel PACS line scan spectra were obtained in Cycle 2@ (P.I. G. Herczeg) in June 
2013 with the goal of searching for mid- to high-J lines of CO, as well as H 2 0 and OH. These 
data cover eight different wavelength regions. The PACS line scan data were reduced from 
level 0 using background normalization scripts in HIPE vll with PACS calibration tree 56. 
Two output spectra were produced, one that was flux calibrated using the central 3x3 spaxel 
region and one that was not. Telescope jitter and slight pointing errors can cause some of 
the source flux in Herschel observations to extend beyond the central spaxel. This effect can 
be corrected by taking the flux ratio of the sum of the central 3x3 spaxels to the single center 
spaxel and comparing it with the ratio of a perfectly pointed observation of a point source, 
then adjusting the source flux accordingly. A comparison between the 3x3 corrected spectra 
and the non-corrected spectra shows a ~10% difference in continuum flux values. Since this 
difference is small, and since the PACS SED spectra were not 3x3 corrected, we chose to use 
the non-3x3 corrected spectra. The continuum flux of the PACS line scan and SED data for 
each region also matched to within ~10%, consistent with the level of systematic error. The 
PACS line scan spectra are shown in Figure [4j 


3. Analysis: Gas Emission Features 
3.1. Emission line identification 

Spectral line identification in both the Spitzer and Herschel spectra was performed both 
by eye and using custom IDL programs. To identify potential lines in the IRS spectrum, we 
first fit the continuum regions of the SH and LH data with a high order polynomial, then 
subtracted the data from the fitted continuum. We then identified candidate emission lines 
as features that span at least two spectral bins and are at least 2cr above the continuum 
(where a is defined as the standard deviation of the continuum subtracted spectra). To 
confirm which of these detections are likely emission lines and not spurious detections or dust 
features, we subtract the continuum dust model (see section 3.1) from the original spectra, 
and repeat the above line detection process. The NIST atomic spectral databas^l, the 


4 Observation IDs: 1342269454 and 1342269455 


5 http : // physics.nist.gov / PhysRefData/ASD/linesJorm.html 
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Splatalogue database for astronomical spectroscopyj, and previous publications pertaining 
to mid- to far-infrared emission features from protoplanetary disks were used to determine 
the species and transition potentially associated with each candidate emission feature. We 
restrict these identifications to spectral features that are at least 2a above the continuum 
and that correspond to molecular or atomic features previously identified in protoplanetary 
disks. Table [1] lists the candidate lines identified in the Spitzer spectrum and their measured 
wavelengths, energies, fluxes, full widths at half maximum, and equivalent widths. 

Emission features in the Herschel PACS and SPIRE spectra were determined in similar 
fashion. The Herschel PACS range scan and SPIRE spectra are less sensitive than the Spitzer 
and PACS line spectra, and thus only a few candidate emission features are identified at levels 
> 3 a. Table [2] lists the candidate lines identified in all Herschel spectra and their measured 
wavelengths, energies, fluxes, full widths at half maximum, and equivalent widths. 


3.1.1. Spitzer IRS molecular & atomic line inventory 


The disk around V4046 Sgr shows evidence for many different atomic and molecular 
species (Figure [TJ. [5] and Table [T]). We detect strong emission from [Ne n], as well as emission 


unco ui ±±2 ^W2 ctiici m tiic spec, 

are commonly seen in T Tauri star/disk systems (e.g. 

u± U.JLJLI. 

Pascu 

cci et al 

20 ll 

3 UI LliCOC O 

Salvk et al. 

JCUCO 

2011; 

PontoDDidan et al. 

2010; 

Pascucci et al. 

2009; 

Salvk et al. 

20(3 I* 

laiita et al. 

2o3. 

The 


Spitzer spectrum also shows H I 7-6 emission at 12.4 prn, which likely does not originate 
within the disk but, rather, arises fro m the accretion shock regions at the (bin ary) central 
stars and/or from their stellar corona flPascucci et al.l 120071: iR.igliaco et al.ll20151 ). 


Emission lines from H 2 O at 15.2 a nd 1 7.2 pm, which are o ften used as markers to 


establish that a disk is water-rich (e.g. Pontoppidan et al. 20101) . are not evident in the 


Spitzer spectrum of V4046 Sgr. The spectrum may display H 2 0 features at 15.4 and 23.9 pm, 
as well as possible complexes of H 2 0 features near 25 and 31 pm (see Figured])- The detailed, 
self-consistent modeling of H 2 0 and OH emission required to confirm these identifications 
will be pursued in followup work. V4046 Sgr also shows HCN emission at 14 .03 pm, but 


lacks the C 2 H 2 emis si on lines typically found in youn ger T Tauri s ystem s flSalvk et al.l 12011 


Carr fe Naiita] 120081) . IPontoppidan et al.l (120101) and ISalvk et al.l (120111) find that H 2 0, OH 
and other molecular gas emission lines originating from warm gas in the inner disk are weaker 
in most transitional disk systems, likely because these disks are beginning to form an inner 
disk clearing either due to planet formation or disk photoevaporation. Therefore, the lack 


6 http : // splatalogue.net/ 
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of strong H 2 O and OH emission lines, and the generally weak levels of molecular emission 
in V4046 Sgr, is consistent with the fact that this system is a transitional disk with an inner 
clearing (IRosenfeld et al. 2013). 


3.1.2. Herschel PACS+SPIRE molecular & atomic line inventory 


Emission features in the Herschel PACS SED and SPIRE spectra and line IDs are 
presented in Figures [2] and [3] and Table [2] The Herschel spectra were obtained with the 
goal of searching for emission from [O i] at 63 /mi and [C i] at 370 pm arising from UV 
photodissociation of CO molecules in the outer disk flMamon et ah 1988). We report the 
detection of [O 1 ] emission (Table [2l Figure [3]), but did not detect [C 1 ]. We also detect 
emission from an OH 3/2-3/2 doublet at 119.2 and 119.4 pm (Figure [3]), which likely arises 
from UV photodissociation of H 2 O in the outer disk. One high level CO ro-vibrational 
transition (J=17 —?• 16) is tentatively detected and 3a upper limits for all other ro-vibrational 
transitions of CO between 65 pm and 600 pm range from 0.42 - 4.68 10~ 14 erg s _1 cm 


-2 


The Herschel PACS line scan spectra cover multiple limited spectral regions also spanned 
by PACS range spectroscopy, but achieved higher sensitivity and spectral resolution. These 
observations reveal several emission features that were below the detection limit of the PACS 
SED data (Figure |T] and Table [2]). Specifically, we detect (or tentatively detect) emission 
from four high-level CO ro-vibrational transitions (J= 31 —> 30, 30 —> 29, 18 —>■ 17, and 15 
—04), as well as OH and CH + . The PACS line scan spectra do not cover the spectral range 
in which emission features were detected in the PACS SED spectra, so we cannot use these 
data to confirm the existence of the H 2 0 and CO J=17 —> 16 features. 


4. Discussion: Gas Emission Features 


V4046 Sgr exhibits emission features that allow us to further understand the struc¬ 
ture, content, and chemistry of the circumbinary disk. We can use these emission features 
to investigate disk photoevaporation and photodissociation processes, and to compare the 
V4046Sgr syste m _with other similarly evolved systems such as TW Hya (age ~8 Myr, d 
~ 53 pc; Ducourant et al. 2013, and references therein), and other T Tauri star-disk sys¬ 
tems in nearby associations. Below, we comment on notable features in the Spitzer and 
Herschel spectra of V4046 Sgr and combine the Spitzer and Herschel spectra to obtain a 
full mid to farUR SED (5-670 pm) so as to compare the SED with a dust disk model from 
Rosenfeld et al. ( 20131) . 
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4.1. H 2 


Although H 2 is the most abundant species in circumstcllar disks, the electric quadrupole 
nature of its rotational transitions result in weak emission that is difficult to detect. However, 


pure rotational H o emission has 


2008 


3 stars (i.e. 

Carr & Naiita 

2011; 

Bary et al. 

200l 

Lahuis et al. 


jeeri detected in the mid- I R around some T Tauri and Herbig 


Naiita et al.ll 2010 l; iBitner et al.l 120081: ICarmona et al. 


2003), and even in a few disks around young brown 


dwarfs (IPascucci et al.ll2013T ) . This emission may be due to collisiona l excitation, X-ray /UV 
irradiation of the disk surface, and/or accretion shocks onto the disk (IBitner et al.ll2008 . and 
references therein). 

We detect emission from pure rotational H 2 0-0 S(l) (J=3—>1) at 17.035_pm in the 
Spitzer spectrum (Tablc[l]). Since V4046 Sgr lies far from any molecular clouds (IKastner et al. 
20081 ). the H 2 emission most likely arises from warm gas in the surface layers of the cir- 
cumbinary disk. The line flux of H 2 is rs-/ 1.50x10 14 erg s 1 cm 2 (L H2 =1.04xl0 28 erg s x ). 
This is comparable to the 17 pm H 2 line strengths measured in T Tauri stars studied by 


Bitner et al. 


(120081 ). as well as to the 17pm H 2 line luminosity from the (si milarly isolated) 


TW Hya (L h 2 =4.17x 10 27 erg s x , assuming a distance of 53 pc; Naiita et al. 2Qloh . 


Carmona et ah (120081 1 did not detect H 2 emission towards V4046 Sgr using VISIR on 


the VLT Melipal telescope. They report a three sigma upper limit for 17 pm H 2 emission 
of <1.3xlO -14 erg s _1 cm 2 , which is roughl y equ ivalent to our measured H 2 17 pm flux. 


The non-detection of H 2 by Carmona et al. (2008) may be due to VISIR having a smaller 
slit size, or the difficulty of removing atmospheric effects from ground-based spectra. It is 
also possible that the H 2 emission is variable as a result of ongoing accretion and variable 
X-ray/UV luminosity ( Barv et al. 2008 1. 


Assuming that the H 2 emission is optically thin, the H 2 gas is in local thermal equi¬ 
librium, and the IRS slit covers the entire source (note that we have already corrected for 
the source possibly extending beyond the coverage of the SH slit), we can estimate an upper 
limit for the mass of the emitting H 2 gas as a function temperature from 


M, as = / x 1.76 x 10 


2n 4:nd 2 XF u i 


hcA u ix u (T ) 


Mo, 


( 1 ) 


where F u i is the flux of the H 2 emission line, A is the central wavelength of the emission line, 
d is the distance in parsecs, A u i is the Einstein coefficient of the J = u — l transition /[ and 


x u {T ) is the population at level u for a given excitation temperature T (Thi et al. 


2001 ). 


19981) 


'A u j=4.76x10 10 s 1 for the IB-ortho J = 3-1 transition (Wolniewicz et al 
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Since only the H 2 -ortho 17 /irn emission line i s detected in our spectrum, we multiply by 
a conversion factor / = 1 + {ortho/par a)^ 1 ([Carmona et al . 20081) . where the equilibrium 
ortho/para ratio can be determined from equation (1) in iTakahash i ( 2001 1. 


Figure [6] shows the resulting mass of the emitting H 2 gas for vari ous exc itation temper¬ 
atures, as well as the total gas+dust mass of the disk established by Rosenfeld e t al. (1201311 
(0.094 M 0 ). It is apparent that if the H 2 gas were at temperatures <100 K, the gas mass 
needed to produce the line flux we see in the Spitzer spectrum would exceed that of the 
total gas+dust disk mass. Indeed, from the upper limit on the flux of the H 2 S(0) emission 
feature at 28.2 //in (<1.9 xlO -14 erg s _1 


cm 2 ), we estimate the lower limit for the excita¬ 


tion temperature of the H 2 emitting gas to be ~100 K (see Thi et al. 2001 . eqn. 3). This 


suggests that much of the IR-emitting H 2 gas resides in the surface layers of the disk, where 
the temperature is warmer than the midplane at a given radius due to stellar X-ray and/or 
UV irradiation. 


4.2. [Ne 11 ] and [Ne in] 


The brightest line in the Spitzer IRS spectrum of V4046 Sgr is that of [Ne 11 ] at 12.8 
jum. We measure a [Ne 11 ] emission flux of 8.67 ± 0.13 xlO -14 erg s' 1 cnr 2 . [Ne 11 ] 
emission is often prominent in T Tauri systems and likely traces hot (T ~5000 K) disk gas 
that is being photoevaporated by EUV and X-ray radiation inci dent on the disk’s upper 
atmosphere within a few tens of AU of the central star (Glassg old et al.l 1200411 . Irradiated 


disk models (IMeijerink et ahl 120081: lErcolano et al.l 12008 


Ercolano fe Owen 2010) predict a 


correlation between the strength of the [Ne ill line and the X-r ay luminosity. Recent studies of 


circumstellar disks (e.g. ISacco et al.l 12012c iGiidel et al.M2010ll have shown a weak correlation 


between X-ray luminosity and [Ne 11 ], and blueshifts of [Ne 11 ] consistent with emission due 
to a photoevaporative wind or magnetically driven outflow. 


Sacco et al. (2012) obtained high resolution rnid-IR spectra covering the [Ne 11 ] (12.81 


//m) emission line towards V4046 Sgr (and 31 other YSOs) using the VISIR spectrograph on 


the VLT. Th ey detected V4046 Sgr and measured a [Ne 11 ] line flux of 6.6 ± 0.2 x 10 14 erg 


— 1 —2 

s cm 


Sacco et ahl (20.12) also measured the flux of the [Ne 11 ] 12.81 /im emission from 
the Spitzer spectrum as 7.68 ± 0.22 x 10~ 14 erg s cnr 2 , which is consistent with our flux 
measurement after rescaling for slit loss. The difference between the Spitzer and VLT [Ne 
II] flux measurements may be indicative of time variability of the emission, but could also 
be due to the different slit sizes. 


The VLT [Ne 11 ] emission from V4046 Sgr is slightly blueshifted and likely arises from a 
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photo evaporative wind is escaping the inner gase ous disk (ISacco et al. 


2012). Rapson et al. 


(2015) show that the R~30 AU gap in the disk (Rosenfeld et al. 2013) is filled with (sub)- 
micron dust grains that may provide enough opacity to produce the observed [Ne ii] line. 
Calculation of the 12 pm opa city throughout the V4046 Sgr disk based on dust column 
densities from Rosenfeld e t al. (120131 ) suggests that there is sufficient optical depth at R~ 
20-50 AU to absorb the red-shifted portion of the [Ne III emissio n. Thus, the [Ne II] emission 
likely originates in this region of the inner disk, as Sacco et al.l (2012) suggested. 


Espai lla t et _ah (120071 ) measure the [Ne ii] emission from CS Cha, a ~2 Myr old T Tauri 


star that is similar to V4046 Sgr in that is surrounded by a transitional disk with an inner 
gap out to R~43 AU. V4046 Sgr has a [Ne ii] luminosity of 5.5 x 10 28 erg s -1 (assuming 
D=73 pc) which is roughly half that measured for CS Cha. If the strength of the [Ne Ii] 
emission is due to X-ray and EUV photons impinging on the disk, then this lower [Ne Ii] 
luminosity would be consiste nt with the fact tha.t the X-ray luminosit y of V4046 Sgr is also 
lower than that for CS Cha (iDonati et al.l 201 ll; Fe ig elson et ah 1993 ). However, both disks 
have a [Ne ii] lumino s ity higher than what is typically observed for classical T Tau r i stars 
(lEspaillat et al.l 120071: IPascucci et al.l 120071 and references therein). lEspaillat et al.l (120071) 
attribute this to the high accretion rate of CS Cha (1.2xl0~ 8 yr _1 ), yet V4046 Sgr has a 

much lower accretion rate (5xlO -10 M 0 yr -1 ; Donati et al. 201 ll) and still shows a high [Ne 


ii] luminosity. Therefore, mass accretion rate may not be a good indication of the strength 
of the [Ne ll] line in circumstellar disks. This comparison points out the need for additional 
studies aimed at investigating potential correlations between [Ne ii] line luminosity and 
star/disk system parameters. 

The 12.8 pm [Ne Ii] feature, combined with measurements of emission from [Ne ill], 
can be used as X-ray/EUV radiation held diagnostics. Emission from 15.6 pm [Ne ill] is 
rarely detected in T Tauri syst ems, and upper limits are usually an order of magnitude 


below that of [Ne ll] (e.g., Lahu is et al. 120071) . We report only an upper limit for the 15.6 


pm [Ne ill] line in Table IR as this line is at best tentatively detected and a potential line is 


only apparent in one nod. For V4046 Sgr, we End [Ne m]/[Ne n] < C 
with the ratio determined for TW Hya ([Ne lll]/[Ne ii] ~ 0.045) (T 

.13, which i 

sfaiita et al. 

s consistent 
2010 ), and 

is comparable to other CTTS that lack jets (e.g. 

Lahuis et ah 

200 

7). Modclii 

lg of X-ray 
'le n] < 0.1 

irradiation of protoplanetary disks ( 

Ercolano et ah 

2008 

) also predicts [Ne m]/[] 


for X-ray bright systems. 


Pascucc i et ah (20.14) constrain the EUV irradiation of the V4046 Sgr disk by studying 


radio data that traces free-free emission from ionized gas caused by EUV irradiation. They 
End that the inferred EUV luminosity reaching the disk is much lower for older disks (V4046 
Sgr, TW Hya and MP Mus) than for younger, less evolved disks, and that the EUV flux for 
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V4046 Sgr is not sufficient to produce the [Ne n] luminosities measured by 
and confirmed in this work. Thus, the disk is likely being photoevaporated by ~1 keV X-ray 
photons at a faster rate than if EUV photons were responsible for the photoevaporation. 


Sacco et al. (2012) 


4.3. HoO and OH 


stars (e.g., 


JLl liUlii ii 2 W CU 

Carr & Naiita 

2011 ; 

11CXD UCC11 UCICCICU 11 

Pontoppidan et al. 

LHC kJ 

2010 ; 

Salvk et al. 

CL CHUU1 

2008; 

iu many yuun^ 

Carr & Naiita 


20081 ). In general, this emission is found to arise from regions of the disk with temperatures 
<1000 K, corresponding to a radius of a few AU, where some H 2 0 is photodissociated into 
OH by far-UV photons. V4046 Sgr shows tentative evidence for H 2 0 emission in the Spitzer 
spectrum, but lacks emission at wavelengths that typically trace the presence of H 2 0 and 
OH in circumstellar disks (see section 3.2.1). In particular, the Spitzer spectrum of TW Hya 
Naiita et al.l (|2010l ) shows no emission from H 2 0 in the Spitzer SH spectrum, but displays 
many OH emission lines. 

Emission from the higher energy rotational states ( E upper >7000 K) of OH observed 
in the TW Hya spectrum likely arises from warm (T~1000 K) H 2 0 gas that has been 
photodissociated within the inner 1 AU of the disk. This sug gests th at, e ven though TW 
Hya ex hibits an inner hole (R~4 AU) in its dusty disk flHughes et ahl 120071: iRosenfeld et al. 


2012bj), there is gas present within this hole in the form of dissociated water vapor. We 
cannot rule out the presence of OH emission in the V4046 Sgr IRS spectrum at the level of 
that in the TW Hya IRS spectrum, because the V4046 Sgr spectrum is less sensitive (due 
to shorter exposure time^l). Thus, while the lack of H 2 0 emission in the Spitzer spectrum 
of V4046 Sgr may be due to FUV radiation photodissociating the H 2 0 to form OH, more 
sensitive measurements and modeling of emission features are needed to confirm the presence 
of OH. 


Emission from OH is clearly present in the Herschel spectra of V4046 S gr, w hile o- 


HoO 


>ia 


—y 7p7 emission at 63 . 32 /ini (Figu re EH) is tent a tively detected. iHoward et al. 


( 2013 1. Riviere-Marichalar et al. ( 2012bh . and Keane et al. (2014) detect ortho-H 2 0 8is —» 
7 o7 emission from various T Tauri disks in the Taurus/Auriga and other nearby star forming 
regions, especially from th ose t hat also display [O i] emission. Based on visual inspection 
of the Herschel spectra in Thi et a l. (120101 1. TW Hya may also show o-H 2 0 emission at 
63//1H. though it is not noted in that paper. Emission from 63 //in o-H 2 0 in T Tauri disks 


8 V4046 Sgr was observed with Spitzer for 30 and 14 seconds in the SH and LH mode, respectively, while 
TW Hya was observed, only in SH mode, for 120 and 600 seconds. 
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had been thought to origi nate from the same regions at which the 


Spitzer spectra is present (IPontoppidan et al.ll2009l; iMeiierink et al. 


but recent modeling of o-H 2 0 emitting sources flRiviere-Marichalar et al 


rot HoO gas emission i n 


2009 


Salvk et al.ll20081) . 


2012 b) has shown 


that the 63 pm emission originates at larger radii in the disk. Similarly detailed modeling 
of the emission spectrum of V4046 Sgr is needed to confirm if this is the case for V4046 Sgr; 
if so, it is a direct probe of cooler H 2 0 in the disk. 


Hogerhei i de et all (20111) detect o-H 2 0 li 0 —> loi and p-H 2 0 In —> 0 0 o emission from 
TW Hya at 538 and 269 pm, respectively, with Herschel-HIFI. They attribute this emission 
to a large reservoir of water vapor at r > 50 AU from the central star. We do not detect 
emission from these H 2 0 transitions in our (relatively low sensitivity) SPIRE spectrum. 

We also detect collisionally excited OH 2 n 3 / 2 emission towards V4046 Sgr at 119.23 and 
119.44 pm in the PACS range spectra, along with radiatively excited OH 2 n 3 / 2 emission 
at 84.60 an d 84.42 pm (blended wit h C O) in the PACS l ine spectra (Figure [3] and HI re¬ 
spectively). Fedele et al. ( 20121 . 2013 ) and Wampfler et al. ( 20131. 2Qlof ) detect H 2 0 and OH 
emission from protoplanetary disks around both low and high-mass young stellar objects. 
Emission from both the collisionally excited 2 n 3 / 2 level, as w ell as the far-infrared radiatively 


excited OH 2 Hi/ 2 level, is seen in most systems. Fedele et al. (2013) also detect OH emission 


features towards protoplanetary disks around both low and high mass stars. Since emission 
from OH 2 Hx/ 2 transitions is not detected in the PACS range scan spectra of V4046 Sgr, 
we conclude that collisionally excited OH molecules may dominate emission around more 
evolved sources like V4046 Sgr, whereas far-infrared pumping and collisional excitations are 
both important processes in younger, less evolved systems. 


4.4. [O 1 ] 


Emission from [O 1 ] is well detected at 63.18 pm in the Herschel /PACS spectrum of 
V4046 Sgr (Figure El Table [2]) - [O 1 ] emission arises in protoplanetary disks when UV 
radiation from the central T Tauri star photodissociates CO or OH molecules in the outer disk 


flMamon et al.l Il988l: lAresu et al.l 120121) . Disk gas mass estimates and other disk properties 


can be derived through observations of 12 CO and 13 CO emission (especially when these 
features are optically thin) so it is important to understand how the CO in a disk has 
been affected by such a UV held. Emission from 10 il has been detected around both 


T Tauri and Herbig Ae/Be stars (e.g 

Keane et al. 

20 14 

; Fedele et al. 2013; Howard et al. 

2013; Riviere-Marichalar et al. 2012a 1 

r, Mathews et al. 

2010 

). and is also well detected in 

the Herschel spectrum of TW Hya (T 

hi et al. 

20101). The [O i] luminositv at 63 pm from 


TW Hya is comparable to the emission we measure from V4046 Sgr, yet neither V4046 Sgr 
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nor TW Hya show detectable 145 /irn [O i] emission or 158 n m [C il] emission flThi et ah 
20101 ). Emission from 370.3 /irn [C i] is also not detected in the V4046 Sgr, but a 3<r upper 
limit is reported in Table [21 This may suggest that UV radiation from the central stars in 
both syste ms is actively photodissociating CO and/or OH in the outer layers of the disk 
flAcke et ahl 2005 1). with the excess C and/or H then being bound up in the form of carbon 
or hydrogen-bearing molecules. 


Keane et al. (2014) used Herschel PACS to search for [O i] 63 /irn emission from full 


disks (i.e. systems with no evidence for a disk gap), transition disks, and outflow systems 
in star forming regions. They detect [O i] emission from 21 transitional disks and find that 


the stre ngth of t 


re [Q i] lines is typically ~l/2 that of full disk sources. Of the stars in 


the Keane et al. (2014) sample, V4046 Sgr is most similar to the full disk systems DK Tau, 


DQ Tau, DS Tau, DG Tau, and HK Tau in Taurus in terms of spect ral typ e, multiplicity, 
and accretion rate. Assuming a distance to Taurus of 137 pc ( Torres et ah 2008i ), we find 
that the luminosity of the [O i] line for V4046 Sgr is ~l/4 that of these Taurus star/disk 
systems, consistent_with the notion that transitional disk systems tend to have lower [O i] 
line strengths (IKeane et ahl 2014). 


4.5. HCN 


HCN serves as a tracer of carbon-rich disk regions, where there is an excess of carbon to 
bond with free H and N atoms. Normally, most of the carbon is bound in either CO or CO 2 


molecules, out 111 cuski 
leaving behind excess 

wit 

c E 

n a mgn n 2 u 

\gundez et al. 

conten 

2008). 

seen in T Tauri stars ( 

Pontoppidan et al 

. 2009) 

rate and X-ray luminosity 

(Teske et al. 2 

011 

)• \ 


Emission from mid-IR HCN is commonly 
and is potentially correlated with accretion 
/Ve detect emission from HCN at ~14 /irn in 
the Spitzer spectrum of V4046 Sgr (Figured]). The integrated HCN line flux is com parab le 


to HCN emission from classical T Tauri stars in Taurus (3.2 mJy /irn; Teske et al. 


2011). 


HCN has previously been de tected aro und V4046 Sgr a t submm wavelengths with the IRAM 
30m and APEX telescopes (IKastner et ahl 120081.12 0141) and this submm emission may also 
be correlated with UV/X-ray luminosity (IKastner et al.l 2 0081 ). 


4.6. CO 


CO emission from the outer (R > 30 A U to ~ 370 AU) disk o f V4046 Sgr has previously 
been detected in the submm with the SMA jRodriguez et al. 2010 ) and IRAM ( Kastner et al. 
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2008 

), and was extensively modeled by 

Rosenfcld et al. 

(2013) 


The Herschel PACS data 
extend these CO detections to higher J transitions that probe the warmer inner regions of 
the disk (R ~10-30 AU). The PACS line scan data (Figure HD show emission from liigh-J 
transitions of CO that are not apparent in the less sensitive range scan PACS data (Figured- 
Specifically, we detect emission from CO J= 31—>30 (blended with OH), 18—>17, and 15—>14, 
and tentatively detect emission from CO J=30—>29. Modeling of these Herschel /PACS CO 
emission features, along with the other emission features mentioned above and previous 
studies of CO emission, would help constrain the gas temperature distribution throughout 
the V4046 Sgr disk. 


5. Silicate Dust 


5.1. Dust grain modeling 

Silicate grain spectral emission features, including amorphous silicates, as well as crys¬ 
talline grain species such as forsterite, enstatite, silica , pyroxene and olivine, have been identi¬ 


fied in the mid-IR spectra, of many T Tauri stars (e.g. ICohen fe Wittebornlll985l: iHonda et al. 


2003; 

Sargent et al. 

<D 

O 

O 

CN 

2009a) 


and, in particular, the type and relative concentrations of silicates in these disks. We mod¬ 
eled the silicate dust features in the S pitzer IRS spectrum of V4046 Sgr using custom IDL 
programs developed by iSargent et al . fl2009bi ). In this method, one fits a two-temperature 
model that is the sum of optically thick emission from warm and cold blackbodies (T^ and 
T c , respectively), and optically thin emission from two dust populations corresponding to 
T w and T c , respectively. The model flux is given by 


F„(X) mod = B v { A, T c ) [fi c + ^ acjKji A) 

j 

A, T w ) Q w + 5 >^(a) 


+ 


( 2 ) 


where B„(A, T) is the Planck function, (H^) is the solid angle of the blackbody at tem¬ 
perature T c (Tu,), a c j (a w j) is the mass weight of silicate dust feature j at temperature T c 
(T w ), and fiy(A) is the opacity at wavelength A for dust species j. While a multi-temperature 
model would be more realistic, a two-temperature dust model that characterizes the inner 
and outer disk regi ons ge nerally well describes the thermal IR emission from T Tauri star 
disks (ISargent et al.ll2009a h To determine the best fit model to the V4046 Sgr spectrum, we 
minimize y 2 between the model and the 7.7 - 33 /mi region of the spectrum. It is assumed, for 
simplicity, that all data points are independent (i.e., the covariance between pixels is zero). 
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Our Spitzer data is excessively noisy longward of 33 //rri, so we omit these data from the 
fit. The best fit model (Table [3]), which has 18 free parameters fitting 1490 individual data 
points, is shown in Figure d This model spectrum assumes a central source of blackbody 
temperature 4200 K, and indicates that the temperature of dust detected by Spitzer lies in 
the range ~120-340 K. 


5.2. Dust Composition 


Table [3] lists the main results of the model fitting just described, in the form of the 
percent by mass of warm and cool dust grain species responsible for the mid-IR emitting 
region of the disk around V4046 Sgr. Note that this model does not include carbon or other 
featureless dust grains that contribute to the continuum emission. The modeling indicates 
that the silicate dust in the inner disk is primarily (86% by mass) composed of large (~5 
/irn radius) amorphous pyroxene and olivine grains at both warm and cool temperatures. 
There are also clear signatures of warm and cool forsterite and silica in crystalline form; 
these grain types constitute the other 14% of the total silicate dust composition, according 
to the model fitting. The modeling indicates that contributions to the mid-IR emission from 
small grains below the Rayleigh limit (27ra/A -C 1, where a is the radius of the dust grain) 
or other non-crystalline components are negligible. 


The foregoing results allow us to compare the composition of the dust in the circumbi- 
nary disk orbiting V4046 Sgr with other protoplanetary disks and other transitional disk 
objects of similar age, such as TW Hya and Hen 3-600 (TWA 3). Modeling of the Spitzer 
IRS spectra of TW Hya implies a large silicate grain fraction of ~34.4%, with only 2.5% 
crystallinity, while modeling of Hen 3-600A results in a similar large si l icate grain percent - 
age of 31.9%, but with a high crystallinity of 36.2% ( Sargent et al. 2006 : Honda et al. 2003 1. 
Studies of T Tauri stars in the Taurus/Auriga star forming region and the TW Hydrae as¬ 
sociation have shown that these transitional disk syste ms have negli gible crystallinity and 
only modest large silicate grain fractions (<35%: ISargent et al. 2006), yet the disk around 
V4046 Sgr has a comparatively high crystallinity fraction and is otherwise dominated by 
large silicate grains. 


Sargent et ah (I2009al ) use methods identical to those described above to model the 


Spitzer spectrum y)f 64 protopla netary disks in the Taurus/Auriga star forming region (age 
~l-2 Myr; Kenyon & Hartmann 1995.). which is considerably younger than V4046 Sgr. They 
find median mass fractions for warm and cool crystalline silicate dust of 11% and 15%, 
respectively, whereas V4046 Sgr displays 3% and 1 4% crys tallinity, r espe ctively. However, 
the transitional and pre-transitional sources in the ISargent et al . (I2009af ) sample all have 
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negligible crystallinity fractions. The fraction of large, cool dust grains inferred for the 
V4046 Sgr disk is also very high- larger than 93% of the sources in Taurus/Auriga. This 
suggests that grain g rowt h has o ccurred in the disk, which is consistent with the advanced 
age of the system. Sargent et al, (2009aj) also find that disks orbiting multiple star systems 
tend to have larger warm large grain fractions than single star systems. This may be due 
to dynamical interactions between the stars and the disk that trigger grain growth in the 
inner disk and/or to photoevaporation of small grains. Since V4046 Sgr is a multiple star 
system with a large mass fraction of cool large grains rather than warm large grains, it is 
possible that a planetary companion in the inner disk is clearing away smaller particles and 
enhancing the large grain population at larger (cooler) radii. 

A variety of mechanisms h ave bee n suggeste d as to how crystalline silicates form in 
protoplanetary disks (see, e.g., Sargent et al. 2006). One of the leading mechanisms involves 
amorphous silicates being thermally annealed into crystalline silicates as a result of heating 
induced by shock fronts within the disk. These shock fronts can result from star-disk inter¬ 
actions, if the disk is massi ve enough, or from local gravitational instabilities associated with 
planet formation. Rapson et al. (1201 51) present direct-imaging evidence for dust segregation 
by size, as well as dust ring structure, in the inner ~40 AU of the disk around V4046 Sgr, and 
suggest that planet formation has occurred or is ongoing in the disk. This planet formation 
activity, if present, would help explain the combination of moderate crystallinity and large 
dust grains. Grain growth and planet formation may have occurred recently in the disk, thus 
creating a larger than typical crystallinity fraction for a disk with an inner gap and advanced 
age. Crystalline silicates at or interior to these planet formation zones would likely accrete 
onto either the central binary star or protoplanet (s). Since the Spitzer spectra are tracing 
the inner few AU of the disk (see section 4.8) we may be seeing the crystalline silicates in 
the innermost portion of the disk that have not yet accreted onto the central stars. 

Another possible scenario for the presence of crystalline silicates in the disk, and the 
varying crystallinity fraction amongst youn g stars, is thermal annealing of surface grains 
due to intense pre-MS stellar outbursts (Abraham et al.1 20091). Such outbursts are typically 
associated with episodes of dramatically enhanced pre-MS accretion rate, however, and the 
(highly evolved, low-accretion-rate) V4046 Sgr system is unlikely to have experienced such an 
episode. Furthermore, our modeling shows that the warm dust is currently at a temperature 
of ~350 K, far below the ~700 K required to thermally anneal dust, so if crystalline silicates 
formed via irradiation during an accretion-related outburst, they have since cooled and 
migrated outward in the disk. 
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5.3. Amending the Rosenfeld et al. (2013 ) model 


Modeling of the entire V4046 Sgr disk was conducted by Rosenfeld et al. (20131) based 
on their SMA CO data, as well as on Spitzer IRS spectra and archival photometric data of 
V4046 Sgr. Using a 3D radiative transfer code, they developed a three-component model that 
includes an inner “gap” filled with pm-sized dust grains out to ~ 30 AU, a ring of cm/nmi- 
sized dust grains from ~30-45 AU, and an extended halo of CO and small grain emission. 
Their resulting model well reproduces the mid- to far-IR photometry and 1.3 mm 12 CO and 
continuum interferometric observations of V4046 Sgr. The same model fits our Herschel data 
rem arkab ly well, due to their inclusion of pm-sized grains out to ~ 30 AU. However, the 
Rosenfeld et all (120131) model did not reproduce the detailed SED structure at Spitzer IRS 
wavelengths. The mismatch between the model and data in this wavelength range was due 
to the presence of a complex (crystalline plus amorphous) silicate grain mixture in the inner 
disk (see section 4.7) that is not accounted for by the three-component model. 


In Figure[Hl we replace the shorter-wavelcngth (7.7-33 pm) portion of the Ro senfeld et al 


020131) model with our silicate model. This revised model, which better reprod uces t he 
Spitze r spectra, merges seamlessly with the Rosenfeld model at ~33 pm. Using the Rosenfeld et al 
0201311 model, we can estimate the characteristic disk radii probed by our silicate model, as¬ 
suming that our newly calculated dust composition and resulting dust opacities do not sig¬ 
nificantly alter the disk temperature in their models. Figure 0 shows disk temperatur e versus 


radius for both large (nun-sized) and small (micron-sized) grains based on the Rosenfeld et al. 


(2013) model. Our silicate dust modeling shows dust at 340 K and 116 K, which correspond 
to a radius of 0.3 AU and 1.3 AU, respectively, for small (pm-sized) grains. We conclude 
that the Spitzer data likely probe the dust grain populations of the inner disk within these 
approximate radii. 


6. Conclusions 


We have presented an analysis of Spitzer and Herschel spectra of the disk around V4046 
Sgr that elucidates the gas and dust constituents within t he disk. We confir m the measure¬ 
ment of strong [Ne n] at 12.8 pm previously reported by Sacco et ah (2012) and report an 
upper limit for [Ne III] emission. The [Ne IIl]/[Ne ll] ratio is consistent with other CTTS 
systems that are X-ray bright and lack jets. We also detect emission from pure rotational 
H 2 0-0 S(l) (J=3—M) at 17 pm that likely originates in the surface layers of the disk where 
the temperature is warmer than the midplane at a given radius due to stellar UV/X-ray 
irradiation. 




















18 


The X-ray/UV radiation coming from the central binary is likely photodissociating CO, 
H 2 0, and OH, resulting in emission from [O i] and OH. We detect emission from [O i] at 
63 /irn, but do not detect [C i] or [C II] in the mid- to far-infrared. Thus, the excess C 
atoms may be forming HCN and other hydrocarbons throughout the disk, and/or the [O i] 
is arising from photodissociated OH in the disk. The strength of the [O i] emission is lower 
than that for star-disk systems of sim ilar spectral type an d ac cretion propert i es sur rounded 
by continuous disks, as was found by I K eane et al. (2014) and Howar d et ah ([2013}). V4046 


emission 


Sgr exhibits OH 2 n 3 / 2 emission at 84.4, 84.6, and 119 /irn, and possible weak H 2 0 
features in the mid- to far-infrared. This suggests that H 2 0 and collisionally excited OH 
molecules reside in the the outer layers of the V4046 Sgr disk. 

Modeling of the Spitzer spectra reveal that the mid-IR emitting region of the disk 
consists primarily of large (~ 5/irn) silicate dust particles (86% by mass), with crystalline 
forsterite and silica particles making up the other 14% of the mass. Overall, the presence of 
large grains suggests that grain growth may be occurring in the inner ~30 AU gap around 
V4046 Sgr and the abundance of cool dust suggests that dust may be settling towards the 
midplane where temperatures are lower at a given radius. The moderate crystallinity fraction 
is simil ar to that of less ev olved protoplanetary disks in the Taurus/Auriga star forming 
region (ISargent et al.ll2009al) . and may result from heating and thermal anneal ing of micron- 


sized dust that is associated with planet formation activity in the inner disk (Rapson et al. 


2015). Combining our si 


i cate dust model with the mid-IR to sub-nun three-component 


model of Rosenfeld et ah (120131 ). we find that the dust emission seen in the Spitzer spectrum 


likely originates from regions of the disk interior to 1.3 AU. These results together suggest 
that grain growth might, in part, explain the transitional appearance of this disk. 
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Fig. 1.— Spitzer IRS spectrum of V4046 Sgr with emission features labeled. Top: 5.5-19.5 
/irn spectrum with an inset showing the H 2 0-0 S(l) (J=3—Fl) line at 17.035 /jrn. Bottom: 
19.5-35 um spectrum with emission features labeled. 
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Fig. 2.— Top: Herschel PACS (black) and SPIRE (blue) spectra of V4046 Sgr. The data 
have been rebinned by a factor of 3. 

Table 1. Emission Lines in Spitzer IRS Spectrum of V4046 Sgr 


Species 

Wavelength a 

(pm) 

K 

Flux 

(10 -14 erg s _1 cm -2 ) 

FWHM 

(10 _2 /zm) 

EW 

(10 _3 /xm) 

H i 7-6 

12.37 

154854 

0.61 ±0.17 

0.32 ±0.09 

0.65 ±0.18 

Ne ii 

12.81 

1124 

8.67 ±0.13 

1.51 ±0.02 

10.00 ±0.15 

HCN 

14.03 


4.90 ±0.38 

8.00 ±0.61 

5.90 ±0.45 

co 2 

14.97 


4.31 ±0.25 

2.77 ±0.16 

5.06 ±0.30 

H 2 0 16 -15* 

15.41 

3602 

1.49 ±0.22 

1.80 ±0.27 

1.76 ±0.26 

Ne ill* 

15.55 

927 

< 1.11 



H 2 S(l) 

17.03 

1015 

1.50 ±0.24 

2.00 ±0.33 

1.62 ±0.27 

H 2 0 9 - 8* 

23.86 

1397 

1.79 ±0.51 

2.74 ±0.79 

1.80 ±0.52 


* Tentative detection. 
a Observed wavelength. 
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Fig. 3.— [O i] (top left), o-H 2 0 (top right), OH 2 n 3 / 2 J= | —> f (bottom left) and CO 
J=17—>-16 (bottom right) emission features detected (or tentatively detected) in the unbinned 
V4046 Sgr Herschel PACS SED spectrum with Gaussian fits overlaid (red). 











































































Table 2. 


Emission Lines in all Herschel Spectra of V4046 Sgr 


Species 

Wavelength 21 

e m 

Flux 

FWHM 

EW 


(^m) 

K 

(10 —14 erg s _1 cm -2 ) 

(10 _2 /im) 

(10 -3 /zm) 


PACS SED 

[0 i] 3 Pi ^ 3 P 2 

63.18 

228 

1.71 ±0.43 

1.30 ±0.33 

6.72±1.71 

o-H 2 0 8 - 7* 

63.31 

1071 

1.09 ±0.27 

1.14 ±0.28 

4.181 ±1.03 

OH 2 n 3/2 j=§ - I 

119.25 

121 

1.11 ±0.13 

9.05 ±0.91 

12.74 ±1.28 

oh 2 n 3/2 J=f - | 

119.43 

121 

1.10 ±0.12 

10.92 ±1.20 

12.65 ±1.40 

CO J=17-16* 

153.28 

846 

1.38±0.25 

13.23 ±2.36 

24.05±4.29 

PACS Line-scan 

OH 2 n 3/2 J=| - § ± CO J=31-30 

84.42 


1.03 ±0.05 

3.30±0.15 

6.49 ±0.29 

OH 2 n 3/2 j=| - I 

84.60 

291 

1.69 ±0.05 

3.29 ±0.09 

11.00 ±0.30 

CO J=30-29* 

87.18 

2565 

0.72 ±0.06 

3.23 ±0.26 

5.21 ±0.41 

CO J=18-17 

144.80 

945 

0.52 ±0.22 

8.78 ±0.38 

8.02 ±0.34 

CO J=15-14 

173.64 

664 

0.84 ±0.02 

10.98 ±0.28 

17.12 ±0.44 

o-H 2 0 2 1;2 - l 0 ,i + CH+ 

179.53 

114 

0.23 ±0.03 

6.72 ±0.99 

5.89 ±0.87 

SPIRE 

C I 

370.30 

63 

< 0.64 




Observed wavelength. 
Tentative detection. 
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Fig. 5.— H 2 0 emission features tentatively detected in the Spitzer IRS spectrum of V4046 
Sgr with Gaussian fits overlaid 


Table 3. Model Dust Mass Percentages in the V4046 Sgr Disk 


Dust type 

% by mass 

Cool Large Amorphous Pyroxene 

20.01 ± 5.79% 

Cool Large Amorphous Olivine 

65.93 ± 8.44% 

Cool Forsterite 

6.96 ± 3.24% 

Cool Silica 

7.09 ± 3.54% 

Warm Large Amorphous Olivine 

96.89 ± 15.55% 

Warm Forsterite 

3.11 ± 4.33% 

Cool Crystalline Silicates 

14.05 ± 5.83% 

Warm Crystalline Silicates 

3.11 ± 7.99% 

Total Large Silicates 

85.96 ± 10.20% 
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Fig. 6.— H 2 gas mass inferred from the strength of the 17 fi m S(l) emission line as a 
function of temperature (solid curve). The horizontal dashed lines represent the masses of 
Jupiter, canonical solar nebula, and the total gas+dust mass of the disk around V4046 Sgr 
( Rosenfeld et ah 2013 1. 
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Fig. 7.— Spitzer IRS spectra (black) with best fit silicate model overlaid (green). The cool 
dust components of the model are shown in red, and the warm dust components are shown 
in blue. 
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Fig. 8.— Syitze r and He rschel s p ectra (black) with our silicate dust model (blue) from 
7-33 urn and th e Rosenfeld et al. ( 20131) model (red solid) from 33-600 fim overlaid. The 


Rosenfeld et ah ( 2013 ) model from 0.25-33 /irn is shown as the red d otte d line. Orange 


aster isks are photometric data points fro m the litera t ure (B & V from iHetg et al.l (120001 ). I 
( 2010 b J,H & K from Cutri et al. ( 2003), AKAR.I 65 fe 90 fx m data from 


from 


Messina et a. 


Yamamura et al.l (12010), WISE bands 1-4 data from lCutri fe et al.l (120131) . and 350,450, 800 


& 1100 /iru data from Je ns en et al. (Il996l ). 
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Fig. 9.— Temperature at a given disk radius in the Rosenfeld et ah ( 2013 ) model for large 
grains (R >5 /jrn; black solid line) and small grains (R ~5 /im; blue solid line). The horizontal 
dashed lines at 340 K and 116 K represent the two temperatures of the micron-sized silicate 
dust grains in our model. 








